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Abstract 

Understanding  the  growth  and  tunction  ot  plant  roots  in  cold  climates  is 
critical,  but  examination  of  root  systems  to  elucidate  their  link  to  production 
is  expensive  and  labor-intensive.  Limited  information  is  available  on  the 
root  growth  and  functions  of  plants  grown  in  northern  climates.  Our 
objective  is  to  present  an  overview  of  the  influence  of  physicol  edaphic 
factors  on  plant  root  systems  with  special  emphasis  on  models  that  are 
available  for  simulating  root  growth.  This  report  summarizes  the  impact  of 
the  soil  physical  environment  (soil  water,  soil  temperature,  soil  air,  physi¬ 
cal  impedance,  and  spatial  variability)  on  root  uptake  and  growth.  Roots 
grow  because  new  cells  are  formed  in  the  meristematic  tissue  near  the  root 
tip,  and  these  newly  formed  cells  increase  in  volume,  pushing  the  root  tip 
forward  if  growth  conditions  are  satisfactory.  Rapid  elongation  of  primary 
roots,  combined  with  well-developed  secondary  roots,  allows  the  plants  to 
exploit  moisture  and  nutrients  from  a  greater  soil  volume.  Root  and  shoots 
are  strongly  interdependent.  The  roots  receive  photosynthates  and  growth 
hormones  from  shoots  and  in  return  furnish  water  and  nutrients  to  the 
shoots.  Several  roof  growth  models  have  been  developed  during  the  past 
decade;  however,  none  addresses  the  problems  associated  with  cold 
regions.  The  models  reported  in  the  literature  can  be  classified  as  1 )  simple 
models,  2)  carbon  partition  models,  3)  growing  degree  day-based  models, 
4)  soil  parameter-based  models,  and  5)  arctic  plant  growth  models. 


For  conversion  of  SI  units  to  non-SI  units  of  measurement  consult  Standard 
Practice  for  Use  of  the  International  System  of  Units  (SI),  ASTM  Standard  E380- 
93,  published  by  the  American  Society  forTesting  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 
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INTRODUCTION 

Plant  roots  grow  because  new  cells  form  in  the 
meristematic  tissues  near  the  root  tip  and  the 
newly  formed  tissues  expand  in  volume,  pushing 
the  root  tip  forward.  The  amount  of  water  and  the 
nature  of  nutrient  uptake  from  the  root  zone  are 
established  by  the  interactions  of  various  physi¬ 
cal,  chemical,  and  biological  processes.  After 
many  years  of  scientific  discoveries,  we  under¬ 
stand  how  water  and  chemicals  enter  the  soil  and 
are  transported  through  the  root  zone.  However, 
our  understanding  of  the  fundamental  aspects  of 
growth  and  proliferation  of  plant  roots  within  the 
cold  soil  environment  is  still  incomplete,  because 
of  the  vast  complexity  of  the  various  interacting 
processes  involved. 

During  the  last  decade,  root  research  has  re¬ 
ceived  its  strongest  emphasis  from  plant  physiol¬ 
ogists,  ecologists,  and  soil  scientists.  Limited  in¬ 
formation  is  available  on  water  and  nutrient 
uptake  by  roots  in  cold  regions.  Billings  (1992)  re¬ 
ported  that  in  northern  climates,  plant  metabo¬ 
lism,  growth,  and  reproduction  are  controlled  by 
the  interactions  of  the  soil  physical  environment 
and  the  genetic  structure  of  cold-adapted  flora.  In 
these  cold  regions,  the  major  constraints  to  root 
growth  are  the  presence  of  continuous  perma¬ 
frost  in  the  subsurface  and  the  short  growing  sea¬ 
son  (Fig.  1),  The  literature  provides  information 
on  the  genetic  structure  of  cold-tolerant  plants; 
however,  limited  information  is  available  on  the 
root  growth  and  development  in  northern  cli¬ 
mates.  Therefore,  this  review  focuses  on  the  soil 
physical  environment  and  its  effect  on  root 
growth  in  cold  climates. 

The  soil  physical  environment  will  be  charac¬ 
terized  by  temperature,  water,  aeration,  mechani¬ 


cal  impedance,  and  spatial  variability.  The  soil 
physical  environment  required  for  the  develop¬ 
ment  of  an  efficient  root  system  must  be  free  from 
stresses  caused  by  mechanical,  water,  oxygen, 
and  cold-temperature  conditions.  Soils,  because 
they  vary  texturally,  structurally,  in  porosity,  and 
pore  size  distribution,  differ  greatly  1)  in  their  re¬ 
sistance  to  root  penetration,  2)  how  they  collect, 
store,  and  release  water  for  root  uptake,  and  3) 
how  they  maintain  aeration  by  gas  exchange  with 
the  atmosphere. 

The  objective  of  this  report  is  to  present  an 
overview  of  plant  root  growth  and  development 
under  soil  physical  environments  with  special 
emphasis  on  models  that  are  available  for  root 
growth  simulation.  The  physical  edaphic  factors, 
as  they  occur  naturally  or  are  created  by  human 
activities,  and  their  effects  on  root  growth  are  pre¬ 
sented,  beginning  with  brief  background  infor¬ 
mation  on  root  origin  and  distribution,  growth, 
development,  and  functions. 

ROOT  ORIGIN  AND  DISTRIBUTION 

Seeds  germinate  and  primary  roots  originate 
from  the  radicle.  Roots  also  originate  from  stem 
nodes  on  stolons,  bulbs,  and  rhizomes  of  grasses, 
or  from  the  hypocotyl  of  certain  dicotyledonous 
plants  (dicots.  Fig.  2).  As  dicots  grow,  the  laterals 
of  the  first  order  emerge  from  the  main  root  fol¬ 
lowed  by  laterals  of  second,  third,  fourth,  and 
even  fifth  order  (Barley  1970).  In  monocotyledon- 
ous  plants  (monocots.  Fig.  2),  the  seminal  (pri¬ 
mary)  roots  originate  from  primordia  in  the  em¬ 
bryo,  and  crown  (nodal,  adventitious,  or  brace) 
roots  arise  from  foliar  nodes  of  the  growing  plant 
(Russell  1977,  Klepper  1991).  The  growth  rate  of 


Figure  1.  Distribution  of  permafrost  in  the  Northern  Hemisphere  (after  Phve  1983). 

Monocots  Dicots 


Figure  2.  Gross  morphology  of  two  root  types;  a  fibrous 
root  system  (left)  and  a  tap  root  system  (right)  in  soil 
(after  Stem  1982). 


the  primary  roots  slows  down  after  the  initiation 
of  nodal  roots.  When  the  formation  of  crown  roots 
is  delayed/  the  plants  rely  entirely  on  the  seminal 
root  system  for  acquisition  of  water  and  nutrients 
during  the  seedling  establishment  stage. 

There  is  a  continuous  spectrum  between  the 
two  extremes  of  fibrous  and  tap  root  systems  (Ren- 
dig  and  Taylor  1989).  Monocot  plants  have  fibrous 
root  systems/  and  dicots  have  some  variant  tap 
root  system. 

Monocots 

Under  favorable  conditions,  monocot  root  ex¬ 
tension  rates  are  rapid.  During  the  vegetative 
phase  of  plant  growth,  root  depth  and  root  length 
density  (m  roots  m"^  soil)  increase  with  time. 
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However,  even  during  seed  sink  initiation,  roots 
do  not  completely  stop  growing.  Root  length 
density  generally  decreases  with  increased  soil 
profile  depth. 

Seminal  and  crown  roots  are  different  morpho¬ 
logically  and  functionally  For  instance,  in  wheat 
and  grasses,  seminal  roots  go  deeper,  are  thinner 
in  size,  and  are  more  branched  compared  to 
crown  roots  (MacKey  1973).  The  crown  root  sys¬ 
tem  is  more  flexible  regarding  development  un¬ 
der  changed  environmental  conditions  (Tennant 
1976). 

Dicots 

The  diameter  of  dicot  roots  decreases  with  in¬ 
crease  in  order  (Rendig  and  Taylor  1989).  The  av¬ 
erage  functioning  life  of  higher-order  roots  may 
be  10  to  20  days.  Higher-order  roots  form  when 
the  environmental  conditions  are  favorable  and 
die  when  water  and  nutrients  are  depleted  in 
their  immediate  vicinity.  When  soil  water  poten¬ 
tial  decreases  below  -0.2  MPa  for  7  to  10  days, 
about  80%  of  the  roots  can  become  senescent. 

Time  and  environmental  conditions  substan¬ 
tially  influence  the  initiation  rates  of  new  roots 
and  death  of  old  roots.  Therefore,  the  root  length 
density  and  root  length  per  plant  fluctuate  during 
the  season.  As  observed  in  monocots,  the  root 
growth  of  dicots  continues  during  the  reproduc¬ 
tive  phase. 

ROOT  GROWTH  AND  DEVELOPMENT 

Roots  grow  because  new  cells  are  formed  in 
the  meristematic  tissue  near  the  root  tip,  and 
these  newly  formed  cells  increase  in  volume, 
pushing  the  root  tip  forward  if  growth  conditions 
are  satisfactory  (Taylor  and  Brar  1991).  Before  the 
cell  starts  expanding,  the  water  potential  inside 
the  cell,  is  equal  to  the  water  potential  out¬ 

side,  \|/Outside^  yinside  =  \j/p  +  \}/o,  where  \|/p  is 

turgor  potential  and  xj/q  is  the  osmotic  potential. 
Processes  responsible  for  cell  expansion  are  de¬ 
scribed  by  Hsiao  and  Bradford  (1983)  as 

1.  Cell-wall  loosening  or  relaxation,  i.e., 

^utside  ^  ^inside. 

2.  Water  flow  into  the  cell,  diluting  its  solutes; 

3.  Water  influx  causing  the  cell  wall  to  expand 
until  x|/^nside  youtside  nearly  equal. 

These  processes  occur  simultaneously  and  con¬ 
tinuously  in  the  plant  root  system. 

Rapid  elongation  of  primary  roots,  combined 
with  well-developed  secondary  roots,  allows  the 
plants  to  exploit  moisture  and  nutrients  from 


greater  soil  volume  (Harris  1977).  Main  axis  root 
elongation  rates  of  forage  legumes  were  signifi¬ 
cantly  correlated  with  seed  weight  (Brar  et  al. 
1990a).  Cultivars  within  the  species  had  signifi¬ 
cant  impact  on  root  elongation  rates.  Similarly, 
root  growth  also  varied  greatly  among  plant  spe¬ 
cies  (Stone  and  Taylor  1983,  McMichael  et  al. 
1985,  Brar  and  Palazzo  1995).  The  rate  of  root 
elongation  may  range  from  0.5  to  3  cm  day"^ 
(May  et  al.  1965)  or  more  (Taylor  and  Ratliff  1969). 
Elongation  rates  of  main  root  axes,  primary  later¬ 
als,  and  secondary  laterals  of  small  grain  cereals 
grown  in  a  relatively  favorable  laboratory  envi¬ 
ronment  were  20, 0.5,  and  0.1  cm/day,  respective¬ 
ly  (Lungley  1973).  The  root  extension  rates  de¬ 
creased  with  increased  plant  age  and  were 
negligible  in  perennial  crops  (Hillel  and  Talpaz 
1976). 

ROOT-SHOOT  RELATIONSHIPS 

Root  and  shoots  are  strongly  interdependent. 
The  roots  receive  photosynthates  and  growth 
hormones  from  shoots  and  in  return  furnish  wa¬ 
ter  and  nutrients  to  shoots  (Taylor  and  Arkin 
1981).  After  seed  germination,  plants  generally 
contribute  the  major  portion  of  the  metabolites 
towards  root  expansion.  For  instance,  the  radicle 
grows  faster  than  the  plumule  of  wheat  after  ger¬ 
mination;  however,  the  trend  is  reversed  at  later 
stages  of  growth.  Rootrshoot  ratios  differ  depend¬ 
ing  on  environmental  conditions  around  the 
plant. 

Relationships  of  root  and  shoot  growth  for  22 
winter  wheat  crops  grown  in  the  United  King¬ 
dom  are  illustrated  in  Fig.  3.  Dry  matter  is  equally 
partitioned  between  roots  and  shoots  in  the  early 
stages  of  wheat  growth,  but  at  anthesis  the  root- 
per-plant  ratio  has  fallen  to  0,1  (Barraclough  et  al. 
1991),  The  total  shoot  dry  matter  of  18  t  ha“^  was 
accompanied  by  a  root  dry  matter  of  1.5  t  ha"k 
Winter  crops  with  longer  growing  seasons  pro¬ 
duced  more  shoot  dry  matter  compared  with 
spring  crops,  and  similar  trends  were  observed 
for  the  root  system.  A  well-developed  root  sys¬ 
tem  was  observed  for  winter  wheat  compared 
with  arable  crops  (Table  1). 

Soil  physical  factors  such  as  temperature  and 
soil  water  potential  affected  the  root:  shoot  ratio 
of  sorghum  seedlings  (Brar  et  al.  1992).  The  ratio 
significantly  increased  with  increases  in  soil  tem¬ 
perature  from  15.9  to  25.2°C  and  declined  at  tem¬ 
peratures  of  30.2  and  35.8°C  (Table  2).  Further¬ 
more,  the  ratio  was  greater  with  -0.1  than  -0.03 
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Figure  3.  The  growth  of  roots  and  shoots  of  high-input 
UK  winter  wheat  crops  in  relation  to  the  number  of 
days  from  sowing.  A  =  anthesis  (after  Barraclough  et 
al  1991). 


MPa  soil  water  potential  at  15.9,  20.5,  and  30.2°C. 
Significant  correlations  (r^  =  0.96)  were  obtained 
for  shoot  dry  matter  with  root  dry  matter  of  tall 
and  fine  fescues  (Brar  and  Palazzo  1995).  Further¬ 
more,  significantly  greater  root-shoot  growth  of 
tall  fescue  was  observed  in  silt  loam  compared 
with  sandy  soil.  Cool-season  grasses  (C3  photo¬ 
synthesis)  had  lower  root:shoot  ratios  compared 
with  warm-season  grasses  (C4  photosynthesis). 

During  the  last  seven  years,  plant  water  rela¬ 
tion  concepts  have  changed,  and  roots  rather 
than  shoots  are  considered  to  regulate  plant 
growth  and  development  in  drying  soil.  Kramer 
(1988)  reported  that  50  years  of  progress  in  our 
understanding  of  plant  water  relations  were  lost 
by  shifting  emphasis  away  from  the  soil  to  the 
shoots  of  the  plant.  Recent  evidence  suggests  that 
root  signals  influence  stomatal  behavior,  leaf  ini¬ 
tiation,  leaf  expansion,  and  other  development 
processes  (Davies  and  Zhang  1991).  Current 
thinking  in  this  area  is  summarized  in  Figure  4. 


Table  1.  Maximum  root  length  density,  MRLD, 
(cm  cm“^)  in  the  topsoil  for  various  field  crops. 


Crop 

MRLD 

Reference 

Winter  wheat 

12.2 

Barraclough  et  al.  (1989) 

Winter  rape 

9.7 

Barraclough  (1989) 

Spring  barley 

4.2 

Welbank  et  al.  (1974) 

Maize  (US) 

4.0 

Mengel  and  Barber  (1974) 

Spring  wheat 

3.3 

Welbank  et  al.  (1974) 

Sugar  beet 

2.6 

Brown  and  Biscoe  (1985) 

Pea 

2.5 

Greenwood  et  al.  (1982) 

Potato  (NL) 

1.9 

Vos  and  Groenwold  (1986) 

Soybeans  (US) 

1.0 

Barber  (1978) 

Broad  beans 

0.8 

Greenwood  et  al.  (1982) 

Table  2.  Root:shoot  ratio  of  sor¬ 
ghum  seedlings  as  affected  by 
soil  temperature  and  soil  water 
matric  potential.  Data  are  means 
of  three  observations  (after  Brar  et 
al.  1992). 

Soil  Root:shoot  ratio 


temperature  Soil  water  potential  (MPa) 
___ra _ -0.03  -0.10 


15.9 

0.34b» 

0.433b 

20.5 

0.39^ 

0.60^ 

25.2 

0.75^ 

0.66^ 

30.2 

0.39'" 

0.44ab 

35.8 

0.36b 

0.06 

LSD  (P  <  0.05) 

0.035 

0.42 

Within  columns. 

means 

not  followed 

by  common  letter  differ  {P  <  0.05)  accord¬ 
ing  to  Fisher's  Least  Significant  Differ¬ 
ence  Test. 

ROOT  FUNCTIONS 

The  major  functions  of  plant  roots  are  water 
transport,  nutrient  uptake,  plant  anchoring,  and 
plant  metabolite  storage.  Roots  absorb  water  and 
nutrients  simultaneously  from  deep  and  shallow 
soil  horizons  and  from  moist  and  partially  dry 
soil  (Rendig  and  Taylor  1989).  Plant  anchoring  is 
important  for  the  shoot  to  emerge  through  crust¬ 
ed  soils  and  for  roots  to  force  a  path  through  a  soil 
matrix  to  keep  the  plant  from  falling.  Storage  ca¬ 
pacity  of  metabolites  is  greater  for  fleshy  roots 
compared  with  grasses.  Roots  are  believed  to  be 
the  primary  source  for  growth  regulators  such  as 
cytokinin,  gibberellin,  abscisic  acid,  and  ethylene. 

Water  uptake 

Soil  water  uptake  is  usually  assumed  to  be 
proportional  to  the  product  of  the  rooting  activ¬ 
ity,  hydraulic  conductivity,  and  potential  energy 
terms  for  any  specific  soil  volume  element  of  imi- 
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figure  4.  Factors  influencing  the  generation  of  chemical  information 
(dotted  lines),  soil  effects  (circles),  plant  physiological  and  developmen¬ 
tal  processes  (rectangles)  in  roots  in  drying  soils  (after  Davies  and 
Zhang  1991). 


form  properties  (Taylor  1983).  A  shallow  root  sys¬ 
tem  with  a  high  root  density  in  the  surface  soil 
can  result  in  rapid  depletion  of  available  water 
and  nutrients,  while  a  deep  root  system  can  make 
available  a  large  volume  of  soil  for  root  extraction 
of  water  and  nutrients  (Chaudhery  and  Sandhu 
1983).  A  deeper  root  system  is  more  efficient  in 
water  uptake  compared  with  a  shallow  root  sys¬ 
tem  (Taylor  and  Klepper  1973). 

In  recent  years,  our  imderstanding  of  the  soil 
water  flow  through  the  soil-plant-air  continuum 


has  been  strengthened  by  analogy  with  electrical 
engineering  concepts  such  as  the  flow  of  electrici¬ 
ty  through  a  network  of  resistances  (Fig.  5).  Axial 
root  resistance  to  water  transport  would  increase 
in  a  deep  root  system  (Meyer  and  Ritchie  1980); 
however,  it  is  an  insignificant  component  com¬ 
pared  with  radial  resistance  (Taylor  and  Klepper 
1971).  Feddes  (1981)  summarized  plant  resistanc¬ 
es  of  various  crops  (Table  3).  These  data  are  empir¬ 
ical  and  can  be  applied  to  conditions  from  which 
they  were  derived.  It  is  evident  from  the  pub- 
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Sink:  Atmosphere  V(t) 


c  ; - Bottom  of  Root  Zone 

•U 


CO 

Q 


Drainage 


Figure  5.  A  root  system  in  layered  soil  viewed  as  a 
resistance  network.  Water  held  in  capacitors  moves 
through  variable  resistance.  Diodes  permit  flow  into  the 
roots  but  not  from  roots  into  the  soil  (after  HUM  1977). 


Figure  6.  Water  flow  in  the  soil-root  continuum.  The 
"valve"  symbols  represent  water  floio  rates  between 
rectangular  storage  compartments  or  to  external  "sink" 
in  the  atmosphere  or  groundwater  (after  Huck  and  HUM 
1983). 


Table  3.  Plant  resistance,  data  of  various  crops  extracted  from  litera¬ 
ture  (after  Feddes  1981). 


Crop 

Type  of 
experiment 

Root  depth 
(cm) 

(10^  cm) 

Source 

Birds  foot  trefoil 

Lab 

60 

11.2 

Gardner  and  Ehlig  (1962) 

Grass 

Field 

26 

10.4 

Rijtema  (1965) 

Potatoes 

Field 

40 

10.3 

Endrodi  and  Rijtema  (1969) 

Red  cabbage 

Field 

42-82 

31.4-36.3 

Feddes  (1971) 

Wheat 

Lab 

45 

9-21 

Yang  and  de  Jong  (1971) 

Wheat 

Lab 

18 

31-58 

Yang  and  de  Jong  (1972) 

Italian  ryegrass 

Lab 

22 

5-17.5 

Hansen  (1974) 

Maize 

Greenhouse 

— 

13 

Reicosky  and  Ritchie  (1976) 

Grain  sorghum 

Field 

100 

11 

Reicosky  and  Ritchie  (1976) 

Douglas  fir 

Field 

80 

25 

Nnyamah  (1977) 

lished  literature  that  measurement  of  resistances 
encountered  in  the  soil-root-plant  system  is  still 
difficult  or  even  impossible  for  many  crops. 

The  pathway  of  water  movement  is  from  bulk 
soil  to  the  rhizosphere,  across  the  rhizosphere 
and  the  soil/root  interface,  across  the  cortex  and 
the  endodermis,  and  then  into  the  lumen  of  the 
xylem.  Water  then  moves  along  the  xylem  to  the 
substomatal  cavities  where  it  undergoes  a  phase 


change  to  vapor  and  diffuses  through  the  stomata 
and  the  leaf/ air  boundary  layer.  A  generalized 
overview  of  macroscopic  water  flow  processes  is 
illustrated  in  Figure  6.  Numerous  water  uptake 
models  are  published  in  the  literature  and  can  be 
classified  as  "macroscopic"  or  "microscopic." 

Macroscopic  models 

Macroscopic  models  assume  the  vertical  tran- 
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sient  flow  of  water  in  a  stable  and  uniform  zone 
of  soil.  The  root  system  is  considered  as  a  diffuse 
sink  for  water  that  permeates  each  depth  of  soil 
layer  uniformly,  though  not  necessarily  at  the 
same  root  length  density  through  the  root  zone. 
Macroscopic  models  are  derived  from  Darcy's 
law,  expressed  in  hydraulic  head  terms 

J^  =  -K^H/dx  (1) 

where  =  flux  of  water  (L^  L“^ 

=  soil  hydraulic  conductivity  (L  t~^) 
dH  =  hydraulic  head  difference  (L) 
dx  =  distance  along  the  flow  path  (L). 

Incorporating  the  soil  water  diffusivity  term, 
D^  =  K^{d^/dQ\  in  eql 

J^  =  ^D^dQ/dx.  (2) 

The  vertical  flow  through  a  thin  layer  of  soil 
where  water  content  changes  with  time  and  dis¬ 
tance  can  be  solved  with  Richards  (1931)  continu¬ 
ity  equation 

dB/dt  =  d/dz  K^dH/dz  +  dK^/dz 

=  d/dz  Dde/dz  +  dK^/dz,  (3) 

Equation  3  is  further  modified  for  water  ex¬ 
traction  by  the  plant  roots 

dO/dt  =  de/dz  (Ky,dH/dz)  +  A(z,t)  (4) 

where  A(z  t)  is  the  root  water  extraction  in  refer¬ 
ence  to  soil  depth,  z,  and  time,  f. 

Several  workers  assumed  that  ^(z,t)  is  a  func¬ 
tion  of  root  activity  None,  however,  has  assumed 
that  root  activity  is  a  function  of  water  potential 
difference  between  plant  and  soil,  distance  be¬ 
tween  the  uniformly  spaced  roots,  and  some 
measure  of  conductivity  in  the  root-soil  system. 
The  root  water  uptake  model  of  Nimah  and 
Hanks  (1973)  can  be  written  as 

7l(z,t)  “  i^roots  ■*"  “  ^(z,t) 

(5) 

where  Hfoots  =  effective  water  potential  in  root  at 
z  =  0 

RRES^z)  =  ^oot  resistance  equal  to  1  +  Rc 

Rc  =  coefficient  to  account  for  longitu¬ 
dinal  resistance  in  xylem 


H(z^t)  “  soil  matric  potential 
S(z^t)  =  soil  osmotic  potential 
RDF (z)  =  proportion  of  total  active  roots  in 
depth  increments 

=  soil  hydraulic  conductivity  at 
depth  z 

Ax  =  distance  between  roots  at  position 
where  H(z^t)  ^rid  S(z^t)  mea¬ 
sured 

Az  =  soil  depth  increment. 

The  major  drawback  of  these  models  is  that  they 
utilize  a  gross  spatial  average  of  matric  and  os¬ 
motic  potentials  and  neglect  the  decrease  in  water 
potential  and  change  in  the  salt  concentration  at 
the  soil/ root  interface  as  well  as  the  rhizosphere. 

Microscopic  models 

Microscopic  models  consider  the  diffusion  of 
water  towards  a  single  root  (Gardner  1960).  The 
models  assume  that  liquid  flow  resistance  in  soil 
is  dependent  on  root  geometry,  rooting  length, 
and  the  hydraulic  conductivity  of  the  soil.  Under 
steady-state  conditions,  the  rate  of  water  uptake 
per  unit  root  length,  from  the  soil  at  a  uniform 
equilibrium  water  content  can  be  estimated  as 

“  “2nR^[Hp  —  Hg]/ln  (^'soil/^root)  (^) 

where  =  soil  hydraulic  conductivity 

Hp  =  matric  head  in  root  epidermis  or  to¬ 
tal  leaf  water  potential 
Hg  =  matric  head  in  soil  surrounding  root 
^soil  =  radius  of  soil  cylinder  surrounding 
root 

rroot  =  ^oot  radius. 

Root  water  uptake  of  a  specific  soil  volume  can 
be  estimated  by  multiplying  the  with  root 
length  density,  Ly.  The  transpiration  rate,  T,  is  as¬ 
sumed  to  be  equal  to  uptake  rate  and  can  be  cal¬ 
culated  as 

7’=i:?ri^vi-  (7) 

i 

Taylor  and  Klepper  (1975)  proposed  the  follow¬ 
ing  equation 

^final  “  ^initial  ”*  (?r)  (^v)  (-^p  ~  ^s)  (^) 

where  Gfinai  and  Ginitiai  are  the  volumetric  water 
contents  at  the  end  and  beginning  of  a  measuring 
period. 

Water  uptake  rates  of  species  differ  even  when 
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Table  4.  Specific  water  uptake  rates  for  some  tropical  and 
temperate  species  (after  Rendig  and  Taylor  1989). 


Species 

Water 
uptake  rate 
(m^  m~^  s"b 

Source 

Cotton 

5  X  10-12  to  5  X  10-11 
<1  X  10-14  to  8  X  10-12 

Bar- Yosef  and  Lambert  (1981) 
Taylor  and  Klepper  (1975) 

Field  pea 

1  X  10-12  to  2  X  10-11 

Hamblin  (1985) 

Lupin 

3x10-12  to  8x10-12 

Hamblin  (1985) 

Corn 

3x10-12  to  2.4  X  10-1“ 

1  X  10-12  to  2.4  X  10-11 

Allmaras  et  al.  (1975) 

Taylor  and  Klepper  (1973) 

Oats 

O 

O 

X 

o 

o 

X 

00 

Ehlers  et  al.  (1980) 

Soybean 

2  X  10-11  to  3  X  10-10 

2  X  10-12  to  3  X  10-10 
2.8  X 10-11  to  2.2  X  10-10 

Allmaras  et  al  (1975) 

Willatt  and  Olsson  (1982) 
Willat  and  Taylor  (1978) 

Wheat 

1  X  10"^^  to  2  X  10“^2 

Greacen  and  Hignett  (1976) 

Spring 

4  X  10“^^  to  lx  10“^2 

Hamblin  (1985) 

Winter 

2  X  10-13  to  2  X  10-12 

Gregory  et  al.  (1978) 

Temperate 

grasses 

about  1  X 

Welbank  et  al.  (1974) 

Figure  7.  Effect  of  plant  age  on  the  hydraulic  conductivity  ofhydro- 
ponically  grozon  and  sand-grozvn  cotton  roots  (after  Brar  etal  1990b). 


they  are  grown  under  identical  conditions  (Table 
4).  For  example,  lupin  and  wheat  were  grown  at 
the  same  location  and  time,  yet  their  maximum 
uptake  rates  were  2  x  10“^  and  1  x  10"^^  m^  m”^  s"^, 
respectively.  Root  uptake  models  often  assume  that 
hydraulic  conductivity  of  the  root  system  is  a  con¬ 
stant  at  all  plant  growth  stages.  Conversely,  Brar  et 
al  (1990b)  reported  that  hydraulic  conductivity  of 
young  cotton  roots  changed  with  plant  age  (Fig.  7). 
Hydraulic  conductivity  decreased  45%  in  nutrient 
solution-grown  roots  and  57%  in  sand-grown  roots 
between  12  days  after  planting  (DAP)  to  the  mean 


of  the  24, 36,  and  48  DAP.  Hydraulic  conductivity 
of  sand-grown  roots  was  4  to  10  times  lower  than 
that  for  roots  grown  hydroponically,  whether 
conductivity  was  compared  at  equal  DAP  or  at 
equal  root  length  (Table  5). 

Nutrient  uptake 

Nutrient  uptake  is  related  to  diversity  in  nutri¬ 
ent  mobility  in  soil  as  well  as  the  surface  charac¬ 
teristics  of  the  roots.  Root  length,  root  radius,  and 
surface  area  are  the  important  characteristics  of 
the  roots  that  determine  the  rate  of  nutrient  up- 
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Table  5.  Total  length  of  cotton 
as  affected  by  plant  age  and 
root  medium;  hydroponic  or 
sand  (after  Brar  et  al.  1990b). 


Root  length 

Plant  Hydroponic  Sand 

age* _ (m)  _  (m) 


12 

3.3+ 

(1.0)*’ 

2.6“ 

(0.7) 

24 

13.4‘> 

(2.4) 

14.4‘> 

(2.9) 

36 

28.8<' 

(1.5) 

17.4+' 

(4.9) 

48 

93.9^ 

(6.2) 

17.1+ 

(2.6) 

*  Days  after  planting, 
t  Means  with  the  same  letter  are 
not  significantly  different  (P  < 

0.05). 

**  Standard  error  of  mean. 

take  (Nye  and  Tinker  1969,  Barber  and  Silberbush 
1984).  The  relationship  between  solution  concen¬ 
tration  {C[ )  of  a  given  nutrient  and  its  uptake  by 
plant  roots  intact  with  the  solution  was  proposed 
by  Bouldin  (1961)  and  modified  by  Nye  (1966) 

/  =  ccci  (9) 

where  }  =  rate  of  nutrient  uptake  (flux)  expressed 
in  terms  of  quantity  per  unit  root  area  per  unit 
time,  and  a=  the  proportionality  term  describing 
the  root  absorbing  power. 

If  the  term  "uptake  per  unit  root  area"  in  the 
above  equation  is  replaced  with  "root  length," 

/  =  lUarCi  (10) 

where  r  =  root  radius.  Nye  and  Tinker  (1969)  com¬ 
bined  the  terms  a  and  r  and  provided  a  new  term 
called  root  demand  coefficient.  For  the  evaluation 
of  total  uptake,  they  incorporated  the  total  root 
length  (L)  in  eq  10 


/  =  2narqL.  (11) 

The  soil  and  plant  factors  affecting  the  root  de¬ 
mand  coefficient  are 

1.  Nutrient  demand  for  plant  growth; 

2.  Operation  of  metabolically  driven  uptake 
and  transport  mechanisms; 

3.  Functioning  of  physical  and  chemical  pro¬ 
cesses  in  the  root  environment. 

Mathematical  models  have  been  developed  for 
nutrient  uptake  by  plant  roots  (Nye  and  Marriott 
1969,  Nye  et  al.  1975,  Claassen  and  Barber  1976, 
Barber  and  Cushman  1981,  Tillotson  and  Wagen- 
et  1982).  These  models  help  in  understanding  the 
complex  soil-rhizosphere-root-plant  system  and 
in  determining  the  significance  of  the  parameters 
involved.  The  basic  transport  equation  used  in 


the  models  to  describe  radial  flux  to  the  root  by 
mass-flow  and  diffusion  as  affected  by  uptake  by 
the  root  is 


dt  r  dr  ^  dr  b 


(12) 


where  Ci  =  ion  concentration  in  solution 

r  “  radial  distance  from  the  root  axis 
ro  =  root  radius 

Dc  =  apparent  diffusion  coefficient 
vq  =  flow  velocity  of  soil  solution  to  root 
b  =  buffer  capacity  =  (3ci/9t)/(9q/3r) 
t  =  uptake  time. 


Boundary  conditions  are  required  to  find  the 
solution  to  eq  12.  The  initial  boimdary  condition 
of  the  equation  is: 


f  =  0,  r>ro,  Ci  =  Cii  (13) 


where  Ca  -  initial  nutrient  concentration  in  the 
soil  solution. 

The  first  boundary  condition  is 


D,b^+ 


at  r  =  ro,  t>0 


(14) 


where  Cmin  =  when  =  net  flux  is  zero,  = 
Michaelis-Menten  constant. 

The  second  boundary  condition,  added  by 
Cushman  to  account  for  competition  between  ad¬ 
jacent  roots,  is 

9c- 

Dj-b +  u^Ci  =  0,  r  =  rj,  f>0  (15) 

dt 

where  q  =  mean  half-distance  between  root  axes. 
In  the  absence  of  competition,  €[  =  Ca  at  r  =  q  and  t 
>0. 

Assumptions  made  in  the  Claassen-Barber 
model  are 

1.  Soil  is  homogeneous  and  isotropic; 

2.  Root  zone  has  constant  soil  water  content 
that  is  maintained  at  field  capacity; 

3.  Plant  or  root  age  has  no  impact  on  nutrient 
flux; 

4.  Nutrient  concentration  has  no  impact  on 
convective  component  of  flux; 

5.  Nutrients  move  to  the  root  by  diffusion  and 
mass-flow  processes;  roots  are  considered 
as  smooth  cylinders  that  absorb  nutrients 
from  the  soil  solution  as  per  Michaelis- 
Menten  kinetics; 

6.  Apparent  diffusion  coefficient  and  buffer 
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capacity  are  independent  of  the  nutrient 
concentration. 

Parameters  required  for  the  Claassen-Barber 
model  along  with  one  additional  parameter  re¬ 
quired  for  the  Cushman  model  are 

1.  Root  morphology 

a)  Initial  root  length 

b)  Rate  of  root  growth 

c)  Average  root  radius 

2.  Root  uptake  kinetics 

a)  Maximum  influx  /max  of  the  nutrient 

b)  External  nutrient  concentration  to  ob¬ 
tain  one-half  /max 

c)  External  nutrient  concentration  Cmin 
where  net  nutrient  uptake  is  nil 

3.  Soil  nutrient  supply 

a)  Rate  of  water  influx 

b)  Concentration  of  the  nutrient  in  the  soil 
solution  at  the  start  of  the  plant  growth 
period 

c)  Apparent  diffusion  coefficient  of  the  nu¬ 
trient 

d)  Differential  soil  buffer  capacity  b  for  the 
nutrient 

e)  Mean  half-distance  between  root  axes 

The  Cushman-Barber  model,  written  in  Pascal 

and  compiled  for  IBM  PC,  is  available  for  educa¬ 
tional  purposes  (Oates  and  Barber  1987). 

ROOT  GROWTH  UNDER  PHYSICAL 
EDAPHIC  CONSTRAINTS 

Soil  water 

The  most  crucial  edaphic  factor  in  a  plant's  life 
is  water,  which  links  it  to  the  soil  via  roots  and 
serves  as  a  vehicle  for  nutrient  transport.  Water 
also  controls  the  exchange  of  gases  and  moder¬ 
ates  soil  temperature  changes  (Clothier  and  Scot- 
ter  1985).  Soil  absorbs,  stores,  and  releases  water 
depending  upon  volume  and  size  distribution  of 
pores  as  determined  by  texture,  structure,  organic 
matter  content,  and  the  depth  of  the  soil.  Plant 
available  water  is  held  in  the  soil  between  -0.01  or 
-0.03  MPa  potential  (in  coarse-textured  and  fine- 
textured  soils,  respectively),  known  as  the  upper 
limit  or  field  capacity,  and  -1.5  MPa  potential, 
called  the  lower  limit  or  permanent  wilting  point. 
Water  availability  may  not  be  limited  to  -1.5  MPa 
potentials;  plant  roots  can  extract  water  at  lower 
potentials  depending  upon  the  plant  type  and  the 
aerial  environment  (Musick  et  al.  1976). 

Root  growth  rates  are  controlled  by  the  pres¬ 
ence  of  continuing  supplies  of  water  to  maintain 
hydrostatic  pressure  in  the  elongating  cells  of  the 


root,  metabolites  for  cell  wall  construction,  and 
growth  hormone  to  loosen  the  bonds  within  the 
cell  wall  constituents  (Lockhart  1965).  Water 
flows  radially  into  the  elongating  root  cells  only 
when  the  cell's  total  water  potential  is  lower  than 
the  combined  osmotic  and  matric  potentials  of 
the  soil  (Rendig  and  Taylor  1989).  Furthermore, 
lowered  soil  water  content  can  shrink  root  diame¬ 
ter  (Cole  and  Alston  1974),  which  reduces  root/ 
soil  contact,  increases  root  senescence,  lowers  soil 
hydraulic  conductivity,  reduces  the  water  poten¬ 
tial  of  the  soil  surrounding  a  specific  root,  and 
decreases  root  hairs,  thereby  affecting  the  water 
and  nutrient  uptake. 

Water  absorption  is  influenced  by  the  geomet¬ 
ric  distribution  of  viable  roots  in  the  soil  profile 
and  the  availability  of  water  at  the  root  surfaces 
(Smucker  and  Aiken  1992).  Absorption  of  soil 
water  is  not  always  linear  with  greater  root 
growth.  Incomplete  root  contact  with  the  soil 
and/or  declining  soil  water  potentials  reduce  the 
water  absorption  efficiencies,  resulting  in  the  pro¬ 
duction  of  excessive  plant  root  surfaces.  Weak 
root  contact  with  the  soil  matrix  results  as  roots 
are  clustered  within  the  macropores  of  soils  and 
smaller  roots  grow  along  soil  aggregates  or  in 
pores  larger  than  the  diameter  of  the  roots  (Lafo- 
lie  et  al.  1991,  van  Noordwijk  et  al.  1992).  Greater 
root/ soil  contact  occurs  as  roots  are  exposed  to  a 
compacted  soil  environment  (Kooistra  et  al. 
1992). 

Lascano  and  van  Bavel  (1984)  calculated  the 
partition  of  the  water  uptake  rate  over  the  root 
zone  as 

=  (16) 

where  Rj  =  rate  of  water  extraction  from  the  root 
zone  or  compartment/ 

'Fsj  =  soil  water  potential  in  that  compart¬ 
ment 

i  =  effective  leaf  water  potential 
RDj  =  relative  root  density  in  the  soil  com¬ 
partment 

PHR  =  specific  plant  resistance. 

The  summation  of  Rj  from  each  soil  compartment 
equals  the  total  water  uptake,  provided  there  is 
no  change  in  water  content  of  the  plant.  In  nature, 
this  is  not  true  because  the  plant  water  content 
changes  diurnally  with  radiation  load  and  chang¬ 
ing  soil  water  content.  Lascano  and  van  Bavel 
also  assumed  that  the  plant  hydraulic  resistances 
are  constant;  conversely,  the  significance  of  plant 
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(18) 


resistance  in  water  transport  was  reported  by  Van 
den  Honert  (1948).  Under  steady-state  condi¬ 
tions,  the  pathway  of  water  movement  through 
the  soil-root-stem-leaves  could  be  expressed  as 
resistances  in  series  (Fig.  5)  using  the  equation 

T=^^s-^r/Rs  =  ^r-WRp  (17) 

where  T  =  transpiration  rate,  4^s,  and  are  the 
water  potentials  in  the  soil  matrix,  at  the  root  sur¬ 
face,  and  in  the  leaves,  respectively;  and  and 
Rp  are  resistances  of  the  soil  and  plant  pathways. 
Rp  is  considered  to  be  dominant  over  Rg  de¬ 
creases  with  increased  transpiration  rate  (Feddes 
1981). 

Steady-state  conditions  rarely  exist  in  the  field. 
Furthermore,  the  plant  root  system  is  dynamic 
(root  senescence  occurs  and  new  roots  emerge), 
root  geometry  is  time-dependent,  and  water  per¬ 
meability  varies  with  position  along  the  root  and 
with  time.  Instead  of  considering  water  flow  to 
single  roots,  a  more  suitable  approach  might  be 
the  macroscopic  one,  in  which  a  sink  term  repre¬ 
senting  water  extraction  by  a  homogeneous  and 
isotropic  element  of  the  root  system  (volume  of 
water  per  volume  of  soil  per  imit  of  time)  is  add¬ 
ed  to  the  conservation  mass  equation  (eq  4). 

Variation  in  root  diameter  has  long  been  ob¬ 
served,  which  further  complicates  the  mathemat¬ 
ical  analysis  of  root  fimction.  Models  of  the  trans¬ 
port  of  ions,  water,  oxygen,  and  other  materials 
from  the  soil  to  the  root  interior  assume  the  root 
has  constant  diameter,  and  transfer  coefficients 
remain  constant  with  time. 

The  simulation  of  water  transport  through  a 
root  system  is  complicated  in  cold  regions  be¬ 
cause  of  the  presence  of  permafrost  and  seasonal¬ 
ly  frozen  soils.  Furthermore,  the  soil  freezing 
front  and  ice  formed  in  the  frozen  ground  en¬ 
hance  runoff  and  reduce  groundwater  recharge. 
A  few  workers  (Taylor  and  Luthin  1978,  Guymon 
et  al.  1980,  Hromadka  et  al.  1981)  studied  the 
complex  processes  of  characterizing  simulta¬ 
neous  heat  and  soil  water  transport  in  a  freezing 
soil  without  considering  vegetation.  The  physics 
of  a  frozen  heterogeneous  soil  profile  include  the 
terms  of  the  soil  energy  balance  and  soil  water 
balance.  The  one-dimensional  energy  conserva¬ 
tion  equation  for  potentially  freezing  soils  can  be 
written  as 

d/dz[h{0)dT /dz\  +  p^C^[dv^T  /3z)  +  Sh  =  C^dT /dt 
-p^L^{ddi/dt)-i-Ly{dp^ /dt  +  d^v/dz). 


The  mass  conservation  equation  can  be  written  as 
d/dz[K{e)  (3 V  /  + 1)]  +  d%  /  p^dz  +  Sm 

=  dQ/dt  +  pidQi  /  pldt  (19) 

where  p£  =  density  of  liquid  water 

C£  =  specific  heat  capacity  of  water 
V£  =  downward  liquid  water  flux 
Cg  =  volumetric  heat  capacity  of  soil 
Pi  =  density  of  ice 
Lf  =  latent  heat  of  fusion 
0i  =  volumetric  ice  content 
Pv  =  vapor  density  in  soil  pore  space 
=  downward  water  vapor  flux  through 
soil 

Sh  =  heat  sink/ source 
Sm  =  water  sink/ source 
0  =  volumetric  soil  water  content 
t  =  time 

T  =  temperature 
z  =  soil  depth 

h  =  thermal  conductivity  of  soil. 

Several  other  workers  used  eq  18  and  19  and 
developed  the  following  types  of  water  transport 
models  (Kung  and  Steenhuis  1986). 

1.  Models  based  on  an  analogy  between  the 
mechanism  of  water  transport  in  unsaturated 
soil  and  in  frozen  soil  (Harlan  1971). 

2.  Models  based  on  the  theory  that  the  pressure 
jump  is  due  to  the  curvature  in  the  liquid 
water /ice /air  interfaces  (Bresler  and  Miller 
1975;  Miller  1973,  1977,  1978;  Miller  et  al. 
1975). 

3.  Models  based  on  the  theory  of  irreversible 
thermodynamics  (Kay  and  Groenevelt  1974, 
Groenevelt  and  Kay  1974,  O'Neill  and  Miller 
1982). 

4.  Models  simulating  temperature  distribution 
in  a  partially  frozen  soil  (Miller  et  al.  1984). 

These  models  were  developed  without  consid¬ 
ering  the  plant  root  system.  During  thawing  sea¬ 
son,  soils  without  vegetation  in  cold  regions  are 
susceptible  to  degradation  processes  of  soil  ero¬ 
sion,  nutrient  losses,  and  organic  matter  deple¬ 
tion. 

Soil  air 

Roots  and  soil  organisms  capture  energy  from 
the  oxidation  of  organic  substances  in  a  series  of 
enzyme-catalyzed  reactions.  Plants  need  molecu¬ 
lar  oxygen  to  respire  and  convert  carbohydrates 
to  carbon  dioxide  and  water.  This  is  an  exother- 
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mic  reaction  and  releases  respiratory  energy  that 
is  utilized  for  many  plant  processes.  The  disap¬ 
pearance  of  molecular  oxygen  triggers  a  se¬ 
quence  of  changes  in  the  physico-chemical  prop¬ 
erties  of  the  soil  (Gambrell  and  Patrick  1978, 
Ponnamperuma  1984).  Such  changes  include  ac¬ 
cumulations  of  reduced  metal  ions,  organic  acids, 
and  volatiles  that  are  potentially  harmful  to  plant 
roots  (Drew  and  Lynch  1980).  Such  accumulation 
to  phytotoxic  levels  requires  time.  The  absence  of 
oxygen  alone  is  sufficient  to  profoundly  alter  the 
plant  metabolism  (Drew  1983). 

In  cold  climates,  cool  autumn  or  early  spring 
temperatures  result  in  a  slow  rate  of  depletion  of 
dissolved  soil  oxygen  by  respiration  of  roots  and 
soil  organisms  (Drew  1992).  The  root  apical  mer- 
istem,  the  zone  of  fastest  oxygen  consumption,  is 
particularly  sensitive  to  oxygen  deficiency  and 
the  sequence  of  metabolic  events  leading  to  cell 
death  (Roberts  et  al.  1984a,b).  In  summer,  crops  of 
cold  climates  consumed  oxygen  at  the  rate  of  1.7 
L  m“^  day“^;  furthermore,  oxygen  consumption 
rates  reduce  to  1.5  L  m“^  day“^  in  winter,  and  to 
about  half  in  the  absence  of  plants  (Currie  1970). 

Sub  optimal  oxygen  concentrations  in  the  soil 
air  occur  because  of  interactions  among  soil  prop¬ 
erties  such  as  porosity,  water  content,  tempera¬ 
ture,  surface  water  movement,  and  continuity  of 
air-filled  pores  with  biotic  activity  (Drew  1983, 
Grable  1966).  Jones  et  al.  (1991)  developed  a  root 
growth  simulation  model  and  considered  the  ef¬ 
fect  of  soil  water  content,  bulk  density,  texture, 
and  plant  genotype  on  soil  air.  Increased  soil  wa¬ 
ter  content  or  bulk  density  reduce  the  oxygen  dif¬ 
fusion  rate  and  affect  air-filled  porosity.  The 
diffusion  coefficient  of  oxygen  in  air  is  approx¬ 
imately  10,000  times  greater  (0.23  cm^  s"^)  than 
that  in  water  (0.26  x  lO-^  cm^  s’l)  at  26°C  (Rendig 
and  Taylor  1989). 

A  mathematical  model  for  the  relationship  of 
root  respiration,  soil  oxygen  supply,  and  internal 
oxygen  supply  was  first  developed  by  Luxmoore 
et  al.  (1970).  Measurements  used  to  predict  the  ef¬ 
fect  of  a  reduced  soil  oxygen  supply  on  root 
growth  and  function  are 

1.  Oxygen  diffusion  rate; 

2.  Oxygen  concentration  in  the  soil; 

3.  Air-filled  porosity  percentage; 

4.  Redox  potential; 

5.  Air  permeability. 

Soil  air  that  is  low  in  oxygen  can  contain  high 
concentrations  of  ethylene,  which  is  a  growth  reg¬ 
ulator  rather  than  a  toxin.  Elevated  concentration 
of  ethylene  inhibits  root  elongation,  cambial 


activity,  and  stellar  differentiation  in  roots  (Ren¬ 
dig  and  Taylor  1989).  An  ethylene  concentration 
of  <  0.4  ppm  can  inhibit  nodule  formation  on  le¬ 
gume  roots.  Furthermore,  ethylene  has  some  ad¬ 
vantageous  impact  on  roots  and  stimulates  later¬ 
al  branching  as  well  as  root  hair  initiation 
(Rendig  and  Taylor  1989). 

Soil  temperature 

Temperature  influences  plant  processes  at  the 
cellular  level,  such  as  osmotic  potential,  hydra¬ 
tion  of  ions,  stomatal  activity  and  transpiration, 
Gibbs  free  energy  available  for  work,  membrane 
permeability,  solute  solubilities,  diffusion,  and 
enzymatic  activities  (Voorhees  et  al.  1981).  The  fi¬ 
nal  shape  of  the  root  system  is  determined  from 
root  branching  and  elongation  of  individual  root 
axes.  Several  reviews  indicate  that  most  plant 
species  exhibit  an  optimum  temperature  for  max¬ 
imum  root  elongation  rates  (Nielson  and 
Humphries  1966,  Cooper  1973).  Furthermore,  the 
species  within  a  genus  (Heinrichs  and  Nielsen 
1966)  and  cultivars  within  a  species  Qohnson  and 
Hartman  1919)  varied  in  root  elongation  and 
branching  rates  in  response  to  temperature.  For 
instance,  Brar  et  al.  (1990a)  reported  high  vari¬ 
ability  among  legume  cultivars  in  main  axis  root 
length  at  different  specific  temperatures  (Fig.  8). 

Low  temperature  reduces  water  absorption  by 
increasing  the  water  viscosity  and  decreasing  cell 
membrane  permeability.  Furthermore,  low  tem¬ 
peratures  also  decrease  metabolic  activity  and 
decrease  the  root  growth  (dry  matter  produc¬ 
tion).  The  addition  of  fertilizer,  particularly  phos¬ 
phorus,  may  compensate  to  some  extent  for  the 
reduced  growth  in  cold  soils  (Nielsen  and 
Humphries  1966). 

Snow  cover  has  a  large  influence  on  soil  tem¬ 
perature.  Snow  cover  decreases  heat  conduction 
to  the  soil  surface  (Legget  and  Crawford  1952) 
and  serves  as  a  sink  for  heat  fluxes  at  both  upper 
and  lower  surfaces  of  the  cover  (Granger  et  al. 
1977).  Relationships  between  air  temperature  and 
snow  depth  can  be  used  to  predict  the  minimum 
soil  temperatures  (Fig.  9).  The  pattern  of  seasonal 
warming  of  air  temperature  is  predominantly 
temporal  at  the  ground  surface;  however,  its  ef¬ 
fect  is  temporal  as  well  as  spatial  below  ground 
and  causes  temperature  limitations  on  root 
growth  (Fig.  10).  During  the  growing  season,  as 
the  warming  front  moves  downward,  tempera¬ 
tures  of  the  deeper  soil  layers  become  suitable  for 
root  growth.  The  depth  of  soil  thermally  suitable 
for  appreciable  root  growth  is  a  function  of  both 
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A)  Paraponto  gama  medic 

B)  Mt.  Baker  subterranean  clover 

C)  Chief  crimson  clover 

D)  Kenstar  red  clover 


E)  Vanguard  vetch 

F)  Woodford  vetch 

G)  Austrian  winter  pea 

H)  Maxidor  alfalfa 


I)  Amcio  arrowleaf  clover 

J)  Cahaba  white  vetch 

K)  Hairy  vetch 

L)  Eski  sainfoin 


Figure  8.  Main  axis  root  length  at  ten  days  after  planting  as  a  function  of  tempera¬ 
ture  (after  Brar  et  al.  1990a). 


Minimum  Soil  Temperature  at  3cm  Depth  (°C) 


Figure  9.  Minimum  soil  temperature  at  3-cm 
depth  as  affected  by  minimum  air  tempera¬ 
ture  and  depth  ofsnozo  cover  (after  Shut' gin 
1967). 

soil  temperature  regimes  and  root  temperature 
response  (Kaspar  and  Bland  1992). 

Soil  temperature  usually  limits  root  system  ex¬ 
pansion  and  proliferation.  Root  expansion  is  a 
function  of  temperature-dependent  processes 
such  as  growth  and  development.  Growth  pro¬ 
cesses  include  cell  elongation,  root  length  in¬ 
crease,  and  root  diameter.  Developmental  pro¬ 
cesses  generally  control  growth  duration  of  cells 
(Burstrom  1956,  Beauchamp  and  Lathwell  1966), 


Figure  10.  Seasonal  warming  of  air  temperatures  (5- 
year  average  of  daily  means)  determines  when  temper¬ 
atures  are  favorable  for  plant  growth.  Soil  temperatures 
favorable  for  root  growth  vary  with  time  and  depth  (mod¬ 
eled  thermal  regime)  (after  Kasper  and  Bland  1992). 
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organs,  and  whole  plants  (Beauchamp  and  Lath- 
well  1967).  The  number  of  lateral  roots,  distance 
from  a  root  tip  to  first  visible  lateral  root,  lateral 
root  density,  root  hairs,  and  root  plagiotropism 
(growth  at  angles  from  the  vertical)  are  tempera¬ 
ture-sensitive  (Atkin  et  al.  1973,  Engels  et  al  1992, 
Mozafar  and  Oertli  1992). 

Root  dry  weight  is  the  most  commonly  mea¬ 
sured  parameter  that  follows  a  typical  growth  re¬ 
sponse  curve  (Cooper  1973).  Root  length,  surface 
area,  and  extension  of  roots  into  deeper  soil  layers 
are  important  parameters  that  influence  water  and 
nutrient  uptake.  Present  technology  has  the  capa¬ 
bility  to  measure  these  parameters.  Root  diameter 
is  another  important  parameter  influenced  by 
temperature;  it  substantially  reduces  water  and 
nutrient  uptake  (Barber  et  al.  1989)  as  well  as  the 
amount  of  carbohydrates  required  for  root  growth 
(Fitter  1991). 

Most  plant  growth  models  related  to  water  and 
nutrient  transport  require  soil  temperature  with 
depth  as  input  for  the  model.  Temperature  simula¬ 
tion  models  published  in  the  literature  are 

1.  Mechanistic  models  that  need  detailed  initial 
and  boundary  inputs  (Richtmeyer  and  Mor¬ 
ton  1967,  Wierenga  and  DeWit  1970,  Hanks  et 
al.  1971,  Rosema  1975,  Van  Bavel  and  Hillel 
1976,  Grant  et  al.  1990); 

2.  Functional  models  that  require  weather  sta¬ 
tion  information  and  soil  information  at  one 
depth  (Neild  1971,  Hasfurther  and  Burman 
1974,  Cruse  et  al.  1980,  Gupta  et  al.  1984,  Par- 
ton  1984,  Brar  and  Unger  1994); 

3.  Statistical  models  of  a  regression  type  that 
were  developed  at  a  particular  site  from  a 
limited  data  set  (Meikle  and  Treadway  1979, 
Ghuman  and  Lai  1983,  Langholz  1989,  Dwyer 
et  al.  1990). 

The  timing  of  snowmelt  and  above-freezing 
temperatures  are  the  major  factors  controlling  the 
onset  of  plant  growth.  Root  elongation  begins  as 
soon  as  the  soil  is  thawed  and  elongation  is  possi¬ 
ble,  depending  on  the  depth  of  the  individual  root 
tips. 

Physical  impedance 

Physical  impedance,  sometimes  called  mechan¬ 
ical  impedance  or  excessive  soil  strength,  can  se¬ 
verely  affect  normal  root  growth  patterns.  Such 
impedances  result  from  increased  soil  bulk  densi¬ 
ty,  increased  cohesion  between  soil  particles,  in¬ 
creased  friction  between  soil  particles,  reduction 
in  soil  water  content,  frost-heave  action  of  soil,  and 
presence  of  permafrost  within  the  root  zone.  The 


Figure  11.  Effects  of  pressure  applied  to  the  outside  sur¬ 
face  of  a  flexible  container  on  relative  root  elongation 
and  shape  of  barley  roots  (after  Russell  and  Goss  1974). 


adverse  physical  impedances  can  severely  affect 
the  plant's  ability  to  emerge  from  crusted  soils,  to 
extend  its  root  system  into  unexplored  soil  vol¬ 
umes,  to  transport  photosynthates  from  shoots  to 
roots,  and  to  transport  water  from  roots  to  shoots. 
These  impedances  restrict  below-ground  expansion 
of  root  crops,  e.g.,  turnip,  sugar  beet,  or  radish. 

Under  an  excessive  soil  strength  environment, 
roots  enter  the  soil  volume  where  pore  sizes  are 
larger  than  the  root  tip.  Conversely,  if  the  pore 
sizes  are  too  small  for  the  entry  of  the  main  root 
but  not  for  laterals,  then  laterals  proliferate  and 
produce  a  highly  branched  root  system  (Russell 
and  Goss  1974).  Furthermore,  if  the  pore  sizes  are 
too  small  and  prevent  entry  of  laterals,  a  stunted 
root  system  can  result.  Maximum  pressures  that 
roots  can  exert  on  the  surrounding  soil  matrix 
range  from  900  to  1500  kPa.  Different  species  exert 
different  maximum  root  growth  pressures  (Taylor 
1983).  The  impeded  roots  are  thicker  in  diameter 
and  distorted  in  shape  (Fig.  11). 

Physical  impedances  result  from 

1.  Natural  consolidation  during  soil  forming 
processes; 

2.  Trampling  by  animals,  including  humans; 

3.  Natural  shrinkage  of  soils  upon  drying; 

4.  Soil  response  to  pressures  and  deformations 
imposed  by  vehicle  wheels; 
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5.  Actions  of  overburden  and  water  droplets  on 
water-weakened  aggregates  during  rainfall; 

6.  Frost  heave  actions. 

Physicalimpedance  canbe  measured  with  pene¬ 
trometers,  shear  vanes,  shear  rings,  triaxial  shear 
cells,  and  unconfined  compression  cells.  As  far  as 
we  know,  no  mathematical  model  is  available  to 
predict  soil  yielding  under  load  using  fundamental 
parameters. 

Spatial  variability 

Considerable  research  has  been  reported  in 
the  literature  during  the  past  decade  on  the  spa¬ 
tial  variability  of  soil  physical  properties.  How¬ 
ever,  in  cold  environments  a  computer  literature 
search  has  revealed  limited  published  studies  on 
the  spatial  variability  of  solute  and  temperature, 
and  none  with  respect  to  root  water  uptake. 
Therefore,  this  section  is  limited  to  reporting  sol¬ 
ute  and  temperature-related  variability  studies 
from  cold  environments. 

Spatial  variability  has  a  significant  effect  on 
field-scale  solute  movement  through  the  unsat¬ 
urated  zones.  The  effect  of  spatial  variability  in 
soil  hydraulic  properties  has  been  investigated 
using  either  parametric  models  (Dagan  and  Bres- 
ler  1979,  Bresler  and  Dagan  1981,  Amoozegar- 
Fard  et  al.  1982,  Simmons  1982,  Destouni  and 
Cvetkovic  1989)  or  transfer  function  models  (Jury 
et  al.  1986,  Butters  and  Jury  1989).  Few  workers 
studied  the  temperature  spatial  variability  in  re¬ 
lation  to  soil  water  content  (Vauclin  et  al.  1982, 
Davidoff  et  al.  1986). 

OVERVIEW  OF  ROOT  GROWTH  MODELS 

Several  root  growth  models  were  developed 
during  the  past  decade,  but  none  is  geared  to  the 
problems  associated  with  root  growth  in  cold  re¬ 
gions.  The  models  reported  in  the  literature  are 
classified  as  simple  models,  carbon-partition 
models,  growing  degree  day-based  models,  soil 
parameter-based  models,  and  arctic  plant  growth 
models. 

Simple  models 

Root  distribution  with  depth 

The  model  of  Gerwitz  and  Page  (1974)  is  em¬ 
pirical  in  nature  and  simulates  root  density  as  a 
function  of  soil  depth.  The  model  makes  no 
attempts  to  define  the  relative  age  classes  of  the  in¬ 
dividual  root  axes  in  each  layer  or  the  proportion  of 
those  axes  that  are  first-  or  second-order  laterals. 


Root  depth  loith  time 

Borg  and  Grimes  (1986)  developed  a  model  to 
simulate  actual  root  depth  as  a  function  of  time. 
For  a  given  crop  and  location,  the  model  requires 
an  estimate  of  the  maximum  rooting  depth  to  be 
achieved  and  the  number  of  days  to  reach  matur¬ 
ity. 

Root  depth  with  temperature 

The  model  is  two-dimensional  and  was  devel¬ 
oped  from  data  that  explain  the  effect  of  root  zone 
temperature  on  soybean  root  development  (Stone 
et  al.  1983).  The  model  predicts  the  extension  of  the 
taproot  and  ten  primary  lateral  roots  for  two 
groups  of  soybean  cultivars.  The  simulation  em¬ 
phasizes  the  importance  of  soil  temperature,  but 
ignores  many  other  factors  responsible  for  root 
elongation. 

Carbon  partition  models 

Carbon  partition  models  are  complex  compared 
to  simple  models.  These  models  partition  carbon 
between  the  above-  and  below-ground  biomass. 
The  partitioning  coefficient  is  assumed  to  be  a  con¬ 
stant  fraction  of  the  daily  photosynthetic  produc¬ 
tion.  Following  are  a  few  examples  of  carbon  parti¬ 
tion  models. 

Rootishoot  ratios 

Numerous  studies  have  shown  that  plants  tend 
to  adjust  their  root:shoot  ratio  to  maintain  an  inter¬ 
nal  carbon-nitrogen  balance  favorable  for  growth 
(Davidson  1969,  Reynolds  and  Thornley  1982, 
Skiles  et  al.  1982,  Johnson  1983,  Coughenour  et  al. 
1984,  Fishman  et  al.  1984).  Davidson  (1969)  devel¬ 
oped  a  model  based  on  a  hypothesis  that  a  func¬ 
tional  balance  exists  between  the  size  and  activity 
of  the  shoot  (which  supplies  carbohydrates)  and 
the  size  and  activity  of  the  root  (which  supplies 
water  and  essential  nutrients).  His  model  suggest¬ 
ed  the  partitioning  of  photosynthates. 

Partition  models  were  reviewed  by  Reynolds 
and  Thornley  (1982).  They  developed  a  model 
based  on  the  assumption  that  the  partition  of  new 
growth  is  controlled  by  whole-plant  substrate  con¬ 
centrations,  Their  model  partitioned  the  photo¬ 
synthates  between  root  and  shoot  on  the  basis  of 
the  nitrogenicarbon  ratio  in  the  plant's  labile  pool. 
In  spite  of  good  features  and  the  high  complexity 
of  the  model,  however,  it  has  the  following  draw¬ 
backs: 

1.  The  carboninitrogen  (C:N)  substrate  concen¬ 
tration  ratio  is  always  considered  a  fixed  value 
not  influenced  by  the  external  environment; 
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2.  The  model  is  unstable  in  terms  of  balanced 
growth; 

3,  The  estimation  of  model  parameters  may 
pose  considerable  experimental  difficulties. 

The  drawbacks  of  the  Reynolds  and  Thornley 
model  were  corrected  by  Johnson  (1983).  His 
model  is  based  on  the  partition  of  new  growth 
between  shoot  and  root  in  terms  of  carbon  or 
nitrogen  substrate  levels  and  the  shootiroot  ratio. 
The  model  considers  the  partitioning  parameters 
and  dynamic  responses  to  changing  environment 
and  root-  and  shoot-specific  growth  rates. 

Grass  crop  model 

Johnson  and  Thornley  (1985)  developed  a 
model  that  divides  root  structure  into  four  cate¬ 
gories:  growing,  newly  expanded,  medium-aged, 
and  senescencing  roots. 

SimRoot 

SimRoot  is  a  four-dimensional  (i.e.,  space  and 
time)  data  structure  model  developed  by  Davis 
(1993)  and  used  by  Nielsen  et  al.  (1994)  to  simu¬ 
late  root  growth  based  on  types  of  root  branches, 
branching  angles,  growth  velocities  of  specific 
types  of  branches,  spatial  distribution  of  respira¬ 
tion,  carbon  exudation,  and  biomass  along  root 
axes.  They  emphasized  the  importance  of  consid¬ 
ering  root  architecture  in  relation  to  nutrient  ac¬ 
quisition  efficiency  and  environmental  factors, 
such  as  elevated  CO2  concentration  that  affects 
below-ground  allocation  of  carbon  resources. 

Growing  degree  day-based  models 

WHTROOT  (Porter  et  al.  1986),  a  winter  wheat 
root  growth  model,  was  developed  based  on 
cumulative  thermal  time  with  descriptions  of  the 
extension  and  branching  of  individual  age  classes 
of  seminal  and  crown  root  axes.  The  model  re¬ 
quires  sowing  date  and  average  monthly  mean 
air  temperature  as  input  parameters.  Model  out¬ 
put  provides  maximum  depth  of  penetration  of 
each  age  class  of  root  and  root  length  density  or 
root  weight  in  any  10-cm  layer  of  soil  contributed 
by  main  axes,  first  order,  and  second  order  later¬ 
als  on  any  calendar  date. 

Soil  parameter-based  models 

Root  grozvth  model 

Jones  et  al.  (1991)  developed  a  root  growth 
model  that  simulates  daily  root  growth  of  crops 
grown  in  a  layered  soil.  Model  input  parameters 
are  aluminum  toxicity,  calcium  deficiency,  coarse 


fragments,  qualitative  constraints,  strength,  aera¬ 
tion,  and  temperature.  The  authors  claim  that  the 
model  will  work  for  mineral  soils  of  temperate  re¬ 
gions.  However,  the  model  does  not  account  for 
soil  situations  such  as  variable-charge  ion  ex¬ 
change  capacity,  wetness-dependent  bulk  density, 
soil  fertility  effect  on  root  growth,  freeze-thaw  im¬ 
pact  on  root  growth,  and  root  water  uptake  in  tem¬ 
perate  climates, 

ROOTMAP 

ROOTMAP  is  a  three-dimensional  model  to  de¬ 
scribe  the  growth  and  structure  of  fibrous  root  sys¬ 
tems  (Diggle  1988).  Root  growth  is  described  in 
terms  of  growing  time,  number  of  axes,  initiation 
time  of  axes,  growth  rate  and  branching  charac¬ 
teristics  of  the  roots,  and  characteristics  governing 
the  direction  of  root  growth.  A  root  system  can  be 
regarded  as  the  result  of  the  accumulated  effects  of 
growth  and  branching  responses  of  individual 
root  tips,  the  root  tissues  behind  these  root  tips,  lo¬ 
cal  soil  conditions,  and  overall  state  of  the  plant. 
The  drawback  of  this  model  is  that  it  cannot  simu¬ 
late  the  entire  root  system  of  a  mature  plant. 

Bengough  et  al.  (1992)  modified  the  ROOTMAP 
model  to  allow  the  numerical  output  of  the  data  on 
root  intersections  with  horizontal  and  vertical 
planes, 

Arctic  plant  growth  models 

Numerous  whole-plant  growth  models  for  vas¬ 
cular  and  nonvascular  plants  have  been  developed 
for  arctic  ecosystems.  Reynolds  and  Leadley  (1992) 
reviewed  21  models  based  on  a  minimal  set  of  pro¬ 
cesses  and  variables  that  are  critical  to  predict 
plant  growth  under  changing  climatic  conditions. 
Most  of  these  models  either  ignored  the  roots  or 
provided  limited  information  on  root  growth  sim¬ 
ulation.  Reynolds  and  Leadley  found  greater  vari¬ 
ability  in  whole-plant  growth  simulations.  They 
classified  these  models  into  four  general  categories 
based  on  degree  of  mechanism  and  comprehen¬ 
siveness:  1)  budget,  2)  flux,  3)  semimechanistic, 
and  4)  mechanistic  (Fig.  12). 

Budget  models  are  empirical  in  nature  and  ig¬ 
nore  any  kind  of  mechanism  or  comprehensive¬ 
ness.  These  models  describe  plant  growth  as  the 
difference  between  input  and  output  parameters. 
In  addition,  budget  models  treat  plants  as  "'black 
boxes"  where  plant  biomass  is  lumped  into  a  sin¬ 
gle  box.  Semimechanistic  models  contain  low  to 
medium  comprehensiveness.  The  plant  growth  is 
described  with  a  minimum  number  of  compo¬ 
nents,  and  heavy  emphasis  is  placed  on  phenologi- 
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Figure  12.  Types  of  arctic  plant  growth  models  based 
on  relative  degree  of  comprehensiveness  and  mecha¬ 
nism.  Compr^ensiveness  refers  to  the  number  of  pro¬ 
cesses  considered  (e.g.,  gas  exchange,  growth,  and 
death),  whereas  mechanism  refers  to  the  level  of  detail  at 
which  each  of  these  processes  is  treated  (after  Reynolds 
and  Leadley  1992). 

cal  modeling.  Mechanistic  models  are  based  on 
the  high  mechanism  and  comprehensiveness  that 
provide  very  detailed  information.  Flux  models 
are  a  special  kind  of  semimechanistic  model.  The 
models  contain  a  high  degree  of  mechanism  and 
comprehensiveness  in  terms  of  gas  exchange  per 
se,  but  tend  to  have  low  overall  comprehensiveness 
in  terms  of  other  plant  growth  processes. 

It  is  evident  from  the  literature  reviewed  that 
whole  plant  models  combine  mechanistic  princi¬ 
ples  with  empirical  observations  to  predict  mass 
as  a  function  of  various  edaphic  factors.  Root  up¬ 
take  is  usually  treated  in  a  highly  simplified  sub¬ 
model,  the  root  system  acting  as  a  zero-sink  for 
nutrients,  and  the  uptake  is  controlled  by  soil  wa¬ 
ter  potential  and  transpiration  rate  or  by  diffu¬ 
sion  flux  rate. 


CONCLUSIONS 

Our  understanding  of  the  fundamental  as¬ 
pects  of  growth  and  proliferation  of  plant  roots 
within  the  temperate  soil  environment  is  incom¬ 
plete.  During  the  last  seven  or  eight  years,  con¬ 
cepts  of  plant-water  relations  in  dry  soil  have 
changed  and  consider  the  roots  rather  than  the 
shoots  to  regulate  plant  growth  and  develop¬ 
ment.  Recent  evidence  suggests  that  root  signals 


influence  stomatal  behavior  and  regulate  leaf  initi¬ 
ation,  leaf  expansion,  and  other  development  pro¬ 
cesses. 
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